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Abstract 

The standard electroweak model is extended by means of a second Brout-Englert-Higgs- 
doublet. The symmetry breaking potential is chosen in such a way that (i) the Lagrangian 
possesses a custodial symmetry, (ii) a stationary, axially symmetric ansatz of the bosonic 
fields consistently reduces the Euler-Lagrange equations to a set of differential equations. 
The potential involves, in particular, a direct interaction between the two doublets. Sta- 
tionary, axially-symmetric solutions of the classical equations are constructed. Some of 
them can be assimilated to embedded Nielsen-Olesen strings. From these solutions there 
are bifurcations and new solutions appear which exhibit the characteristics of the recently 
constructed twisted semilocal strings. A special emphasis is set on "doubly- twisted" so- 
lutions for which the two doublets present different time-dependent phase factors. They 
are regular and have a finite energy which can be lower than the energy of the embedded 
twisted solution. Electric-type solutions, such that the fields oscillate asymptotically far 
from the symmetry-axis, are also reported. 
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I. INTRODUCTION 



Since the discovery of topological string solutions in the so called "Abelian-Higgs- 
model" by Nielsen-Olesen [1], more than thirty years ago, there have been several 
investigations for applying these solutions in high energy physics and cosmology. 
The existence of these solutions is closely connected to the symmetry breaking of 
the local U(l)-symmetry by a Brout-Englert-Higgs mechanism. In the last years it 
has been shown that string-like configurations also appear as classical solutions in 
several other models within several pattern of symmetry breaking. One of the most 
famous generalisation of the string is the "semi- local" string 2|. It appears in models 
involving both local and global symmetry breaking, see [3] for a review. One of 
the most striking properties of these solutions is that they constitute stable vortices 
without being supported by a standard topological argument related to the homotopy 
group of the vacuum manifold Q| . The semi- local string appears naturally in models 
involving N charged complex scalar fields coupled to electromagnetism and presenting 
an SU(N) global symmetry suitably broken by a Brout-Englert-Higgs potential. Along 
with the basic Nielsen-Olesen (NO) strings they possess a cylindrical symmetry, they 
do not depend on the coordinate associated to the symmetry axis (usually chosen as 
the 02-axis) and have a finite energy per unit length of the z coordinate. 

Very recently, the classical equations of this model (in the case N=2) were rein- 
vestigated 0, [f| and new families of solutions were shown to exist. These remarkable 
solutions are called "twisted-semilocal strings". The crucial difference between the 
twisted semilocal and the conventional ones 0, 0] resides in the fact that the two 
scalar fields have different phase factors depending linearly on time as well as on the 
coordinate corresponding to their axis of symmetry. This needs an axially symmetric 
ansatz but, e.g. does not invalidate the cylindrical symmetry of the energy density of 
these solutions. 

In this paper, we consider the semi-local strings (twisted or not) in the framework 
of a model involving two Brout-Englert-Higgs doublets. This is motivated by the 
Minimal Supersymmetric Standard Model (MSSM). Along with [1, 0|, we consider 
the classical equations in the limit where the Weinberg angle 9w is set to 6w = 
7r/2. The most general potential used in the context of MSSM is very involved 
and depends roughly on ten parameters. Here we limit our investigation to the 
subfamily of parameters allowing for a custodial symmetry. This restricts the number 
of parameters to five. We give numerical evidence that, in a specific but rather large 
region of the space of the parameters (involving to a large extent the physical region 
determined by the experimental lower bounds for the Brout-Englert-Higgs bosons) 
solutions exist for which the two doublets present independent twists in the phases 
associated with specific components of them. We call this type of solutions "doubly- 
twisted" , in general they coexist with the embedded simply-twisted solutions [1, @] 
where one of the doublets does not depend on t and z. For some region of the 
parameters, doubly twisted solutions have the lowest energy. The parameters labeling 
the twist can, in particular, be tuned in such a way that the associated Noether 
currents have opposite signs. Another feature of doubly twisted solutions is that, for 
generic values of the potential, the components A and A z of the Maxwell field turn 
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out to be not proportional to each other. This contrasts with the case of one doublet 

The description of the model is presented in Sect. II, where a special emphasis is 
set to the parameters of the potential and their usual relations to the masses of the 
Higgs-particles. The ansatz for the fields and the corresponding equations are given 
is Sect. III. In Sect. IV, we characterize the constraints of the potential that allow 
embedded NO strings. This is useful since it leads to a family of semi-analytical 
solutions (in a sense that they are parametrized in terms of the profile of the NO 
solution only). Several features of the solutions constructed numerically for generic 
values of the potential are discussed in Sect. V. Due to the possibilities allowed by the 
twist with respect to the two doublets, the solutions depend on several independent 
parameters. We present a few plots which hopefully capture a large part of qualititive 
properties of the twisted semi-local string available in the MSSM. Some conclusions 
are given in Sect. VI. 



II. THE MODEL 

The lagrangian that we consider in this paper reads 

(i) 

standard definitions are u 

for the covariant derivative and gauge field strengths : 



Where $( 2 ) denote the two BEH-doublets and the standard definitions are used 



F% = W - d u V; + ge ab XK c (2) 

D^=(d»- l -gT«V«- l -g>A,y (3) 

(the case 9w = 7r/2, i.e. g = will be used later in this paper). 

The most general gauge invariant potential constructed with two BEH-doublets is 
presented namely in [7|, |8( and depends on nine coupling constants. Here we consider 
the family of potentials of the form 

(4) 

depending on five parameters. The direct coupling between the two doublets is 
parametrized by the constant A3 and one of the main characteristic about this po- 
tential resides in the fact that it imposes a symmetry breaking mechanism to each of 
the BEH-doublets. Several classical solutions in the models have been studied in the 
past. Sphericall y sy mmetric sphaleron solutions were studied at length in the case 
A 3 = in 0, 10, 11] and a similar study for A 3 7^ was reported in 12 



The lagrangian (CQ) is of course invariant under SU(2) gauge transformations but 
it further possesses a larger global symmetry under SU(2) x SU(2) x SU(2). In fact, 
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the part of the lagrangian ([I]) involving the scalar fields can be written in terms of 
2x2 matrices defined by 

When written in terms of the matrices Mi and M2, the lagrangian ([ID becomes 
manifestly invariant under the transformation 

VJ = AV^A^ , M{ = AMiE , M2 = AM 2 C (6) 

with A,B,C G SU(2); this is the custodial symmetry. The double symmetry breaking 
mechanism imposed by the potential (T4|) leads to a mass My/ for two of the three 
gauge vector bosons and, namely, to two neutral BEH-particles with masses M h , M H . 
In terms of the parameters of the Lagrangian, these masses are given by [s| 



?\Ai 2 + ^2 2 , M% ih = ] : [A 1 + A 2 ±y/(A 1 -A 2 )* + 4B*] (7) 

with 



At = 2v 2 1 {\ 1 + A 3 ) , A 2 = 2^ 2 (A 2 + A 3 ) , £ = 2A 3 ^ 2 (8) 
For later convenience we also define 



tan/3 =— , p Hh = —— , e p = 4— , p = l,2,3 (9) 
v t M w g 2 



Note that the mass ratios pi )2 used in [10| are related to pjj,h by Ph — m ax{p!,p 2 }, 
Ph = min{pi,p 2 }- F° r physical reasons, we consider only v\ > 0, v 2 > so that 
< (3 < 7r/2. Interestingly, the parameter e 3 can be negative but cannot take 
arbitrary values. The following relations are useful to determine the physical region : 

ei cos 2 p + e 2 sin 2 (3 = - (p 2 H + p 2 h ) - e 3 

1 



e lC os 2 /3-e 2 sin 2 /3 = - ^(p 2 , - p 2 ) 2 - 4e| sin 2 (2/3) - e 3 cos(2/5) (10) 
The physical domain is then determined by the conditions 

pl ~ pl <^<-^d , o<p 2 <pi (ii) 



2 sin(2/3) ~ ~ 2 sin(2/3) 



III. CYLINDRICAL SYMMETRY 



In this paper, we will assume that only the U(l)-subgroup of the gauge group is 
really gauged (that is to say that we assume = V 2 = 0) and look for solutions 
presenting a cylindrical symmetry. For this purpose, we impose a form for the fields 
inspired from [5|, 0]. After some algebra, it turns out that the form 

<P m = h{p) exp^ , 2(1) = / 2 (p) exp mi W^ t+ ^> (12) 



4 



01(2) = £i(p)exp in2 *exp i(<T o i +^) ; 2(2) = ^ 2 (p)exp im ^exp i(CT0 ' +<T32) (13) 

A = uj a (p) , A p = , A (j) = na(p) , /4 3 = w 3 a 3 (p) , (14) 

ni,n 2 and mi,m 2 being integers is consistent with the classical equations. As pointed 
out in [H, HI, solutions exist for < mi < n% and < m 2 < n 2 . Inserting this ansatz 
into the energy functional, it turns out that finite energy configurations are possible 
only for n\ = n 2 = n and a' = a' 3 = 0. In fact we parametrize the U(l) field as close 
as possible to [H|, @] (although without losing generality) in order to crosscheck our 
solutions with the solutions constructed of [5j, |6j in the following two limits : (i) the 
second doublet $( 2 ) becomes trivial, (ii) the parameters o"o,3 appearing in this doublet 
vanish. 

The ansatz above leads to a reduced lagrangian density of the form : 



2/ I \2 2/ /\2 ^ 2/ AH 

^oK) - w 3 (a 3 ) - — n (a) 
P 



p(cos/?) 2 [(/o 2 + U2) 2 + (g[) 2 + W2) 2 + {-^A + <4<% + n2{ \ a? )fl 

p 



2,1 i'2 2,, i2 1 ( m l na ) 2 \ r2 1 / 2 2 1 2 2 1 ^(1 a ) 2 \ 2 

+ (-w„(l-a ) + w 3 (l-a 3 ) + — )/ 2 + (-w a + w 3 a 3 + )# x 

2 ,2 ''." i2 ^_na) 2 N 2l 

P 2 

p(cos/?) 4 [ ei (|/| 2 - l) 2 + e 2 (|^| 2 - tg 2 /3) 2 + e 3 (|/| 2 - 1)(|^| 2 - tg 2 /5)] (15) 



2 2 (m 2 — na) 2 

+ (-(O"0 - ^oOo) + (0-3 - ^3«3) H ^ )# 2 



where the prime means the derivative with respect to the axial variable p. The 
quantity Q defined according to 

Q = 2ttJ dppazUl + gj) = 2vr J dpp(l - a 3 )(/ 2 + g 2 2 ) (16) 

is also of great interest. It determines the vortex worldsheet currents and generalizes 
Eqs.(5.4)-(5.6) of [6]. They are associated with the conserved quantities related to the 
Noether current of the residual U(l) symmetry. In passing we note that the equality 
in (fTB"|) provides a useful test of accuracy for our numerical solutions. 

In order to write the corresponding Euler-Lagrange equations, we find it convenient 
to define K = (ao + a 3 )/2, A = (a — a 3 )/2. It is straightforward to establish the 
Euler-Lagrange equations : 

pC-)' = 2(cosP) 2 {f 2 (a-l) + f 2 2 (a-^)+gl(a-l)+g 2 2 (a-^)} 
pi n n J 

\pf'J = h[-^al + u 2 A + n2{l ~ 2 a? 

+ (cos/?) 2 {2e 1 / 1 (|/| 2 - 1) + ^(M 2 - tg 2 /?)} 



kpffl = /2(-^o 2 (l- a o) 2 + ^ 3 2 (l-a 3 ) 2 + (mi ^ 
P I P 
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+ (cos/3) 2 {2 ei / 2 (|/| 2 - 1) + e 3 / 2 (M 2 - tg 2 ft)} 

/ / \ / | 2 2 2 2 ^ v ^0 

P I P 

+ (cos/3) 2 {2e 2 ^(|^| 2 - tg 2 /3) + e 3 ^(|/| 2 - 1)} 
-{pg'2)' = 92 i-(o"o - ^oao) 2 + (<r 3 - w 3 a 3 ) 2 + ( m2 nQ ) 

p I p 

+ (cos/?) 2 {2e 2 £ 2 (M 2 - tg 2 /?) + e 3 ^ 2 (|/| 2 - 1)} 
P I ^0 ^3 

1/_A/V ON 2 J oAl ^12 , OAUI2 „2/°0 ^3 • 



-(PA')' = (cos/5) 2 2A|/| 2 + 2AM 2 - g 2 (^ - ^) (17) 
P { w w 3 J 

where \f\ 2 = / 2 + /| and |g| 2 = g 2 + gr|- The interest for using A appears on the last 
equation; defining R = (ctq/ljq) , R 3 = (ct 3 /uj 3 ) and 5 = R — R 3 , we see that 5 = 
(corresponding to parallel twists in the two doublets) implies A = by the positivity 
argument. If however solutions exist such that 5^0, they will have A^O; this is a 
new feature with respect to the case with one doublet (1,0). 

In order to obtain regular and finite energy solutions, the above system has to be 
solved with the following boundary conditions respectively at p = and p = oo : 

«o(0) = 4(0) = o(0) = /x(0) = f'M = 9l (0) = g'M = (18) 

a (oo) = a 3 (oo) = / 2 (oo) = # 2 (oo) = , a(ao) = 1, /i(oo) = 1, 01(00) = tan(/3) 

(19) 

Before discussing the solutions for generic values of the potential's parameters, it is 
useful to study a few specific limits of these equations. 



IV. PARTICULAR SOLUTIONS 

A. Nielsen-Olesen strings 

Setting e 2 = e 3 = ft = in the potential and a 3 = a = f 2 = gi = 9\ = 0, the 
system of equations above naturally reduces to the two equations : 

1 1 n 2 

a" - -a' = 2f 2 (a - 1) , f" + -f = f-(l - a) 2 + 2e/(/ 2 - 1) , e = e x (20) 
P P P 

whose solution is the NO string. For later convenience, we denote this particular 
solution by a e , f e . 

B. Embedded Nielsen-Olesen string 

The system of seven equations above admits embedded NO solutions of the form 
a = a € , fi = f e , gi = tan(/?)/ e , a = a 3 = f 2 = g% = (21) 
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provided the following relations between the parameters of the potential hold 

2e = 2ei cos 2 (3 + e 3 sin 2 (3 , 2e = 2e 2 sin 2 (3 + e 3 cos 2 /3 . (22) 

Accordingly, embedded semilocal strings will exist if e > 1/2. For coupling constants 
away from these constraints, we expect the solutions to be deformed progressively 
from the embedded NO-solutions. The linearized equations [ij associated with the 
fields f 2 and g 2 corresponding to an embedded NO-solution lead to a Schroedinger 
equation 

--{pf 2 )'+[- 2 {na-m) 2 -2t{l-f 2 )]f 2 = -uj 2 f 2 , m = m, = m 2 (23) 
P P 

and an identical equation for g 2 with the spectral parameter u 2 replaced by a 2 . Using 
the reasoning of p, \^ , we conclude that superconducting solutions with both fields 
$(i) and $(2) twisted will exist as soon as the potential parameters are such that 
2e > 1. These will exist for 

^ 2 G]0,^ 2 ] , a 2 e]0,oo 2 } (24) 

where u 2 are computed in 0] (see e.g. Tables I, II and Fig. 3 in |6]). As a consequence, 
we expect solutions to exist with arbitrary relative phases associated with the two 
Higgs doublets. This is indeed confirmed by our numerical analysis. See e.g. Fig. 
4 where the energy, the charge Q and the parameters A(0),/ 2 (0) are reported as 
functions of the ratio 0-3/003, for 003 = 0.017. 

It is a natural question to ask whether the bifurcation point corresponding to 
2e = 1 in (|22"1) belongs to the physical region of the parameter of the MSSM. In terms 
of the mass ratios pn,Ph, it turns out that embedded NO-solutions exist if 

p 2 H + p 2 h = 2 + e 3 , |e 3 |/(4 + 5cos(2/?) 2 ) 1 / 2 = |p 2 ,-p 2 | (25) 

If we fix ph, the conditions above define a point 

PH = V / 2+^ > e 3 = (4 + 5 2 c ^fi/ 2 - a ' ^ - Si ^3) (26) 

in the Ph,^3 plane. After some algebra, one can check that, if ph < 1, this point is 
inside the physical region for generic values of (3 and on the curve determining the 
physical region for (3 — 7r/4. Since the experimental results show that the Higgs mass 
is larger than the W-boson mass, no NO-embeded solution is available for realistic 
values of the Higgs mass parameter. 

C. Embedded and mixed semi-local strings 

The twisted semi-local strings of 0, 0] are characterized by the winding number n 
and form n families of solutions labelled by m with < m < n— 1. Around p = 0, the 
two components of the Higgs field behave according to f\ ~ fi P n arid f 2 ~ f 2 P m ■ 
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When the condition (1221) holds, these solutions can also be embedded in the system 
above provided g\ = /itan/3 and g 2 = / 2 tan/3 and 6 = 0, implying A(p) = as 
already mentionned. These embeddings holds for any n and m% = m 2 . However, it 
is likely that solutions can be constructed with mi ^ m 2 and 5 = 0, we will see that 
this is indeed the case and call these solutions 'mixed'. 



D. Simply and doubly twisted solutions 

For generic values of the potential, the system of seven equations above admits 
solutions where all the functions are non trivial and the parameters Rq, R 3 are inde- 
pendant. We will denote these solutions doubly twisted ones. However, as can be 
seen from the equations and the boundary conditions, solutions possibly exist which 
have g 2 = and f 2 7^ (or / 2 = and g 2 7^ 0). In the following chapter, we will 
denote such solutions simply twisted. 



V. NUMERICAL RESULTS 



The system above cannot be solved analytically, we relied on numerical method 
to construct solutions; the numerical routine used is based on the collocation method 
20| . Because both twisted and untwisted solutions obey the same set of bound- 



ary conditions, for instance (1181) . (1191) . some trick has to be implemented in order to 
enforce the numerical routine to produce the twisted strings. For this purpose we sup- 
plement the system by three trivial equations, namely duj^/dp = do^jdp = do^jdp = 
and we take advantage of these supplementary equations to impose by hand values 
for / 2 (0), g 2 (0) and 03(0) as supplementary boundary conditions. If we start from an 
appropriate initial configuration, the values of u 3 ,a , a 3 for which these supplemen- 
tary boundary conditions are fullfilled can be reconstructed numerically. Once such a 
solution is available, the trick allows to "navigate" very efficiently into the manifold 
of simply or doubly twisted solutions. 



A. One doublet 

In order to test the efficiency of our technique, we first reconsidered the one-Higgs 
equations in the case e\ = 2, e 2j3 = (corresponding to (3 = 4 in the notation of [6[]). 
The second doublet then vanishes identically and the equations reduce to the ones of 

@]). The numerical values that we obtain are in full agreement with those of 0. 
Several parameters characterizing this family of solutions are represented on Fig{TJ 



B. Embedded solutions 



Assuming i? 3 = R , we manage to embed the solutions of [6| in the MSSM model 
discussed above i.e. solutions with f\ = gi, f 2 = g 2 . On Fig. 2, we present several 
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FIG. 1: The energy, the charge Q and parameters 03(0), u are plotted as functions of the 
value /2(C)) for the solution with one doublet and ei = 2. 



parameters characterizing these embedded solutions for t\ = e 2 = 2 £3 = and (3 = 
7r/4 and for the solutions corresponding to n = 1, m\ = m2 = 0; n = 2, mi = m2 = 
and n = 2, mi = m 2 = 1. On these plots, the different parameters are represented as 
functions of 10. The plots confirm that the embedded solutions (like the original ones) 
exist only for a finite interval of uo and that the solution bifurcates into an embedded 
NO solution (with f 2 = g 2 = = 03) in the limit uo — > uob- In the limit u — >• 0, a so 
called Skyrme type solution is approached but we did not study this limit in details. 
We just note that the charge Q becomes infinite in this limit, as a consequence of 
the fact that the function a 3 reaches its asymptotic value 03(00) = more and more 
slowly (see 0] for more details). 



C. Mixed solutions 



The profile of a mixed a solution is presented on Fig. [3] corresponding to n = 
2, mi — 0, m 2 = 1, 5 = and /2(C)) = 0.8; a branch of solutions of this type can be 
constructed by varying /2(C)) A few parameters characterizing this branch with respect 
to uo are represented in Fig. 2. For the parameters choosen, it exists for uo < 0.32. In 
the limit uo — > 0.32, we have g' 2 {0) — > and the function g 2 tends uniformly to zero 
while 03, f 2 remain non zero, this feature is specific to mixed solutions. 
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FIG. 2: The energy, the charge Q and parameters 03(0), /i(0), /^(O), etc... are plotted as 
functions of u for the embedded solutions corresponding to a few values of n,mi,m2 ■ 

D. Doubly twisted solutions 

We first investigate the effects of a doubly twisted solution of the NO-string type 
corresponding to e\ = e% = 2, €3 = 0, (3 = 7r/4 and assuming u 3 = 0.017 , Rq = 1. 
The effect of the parameter R 3 on the solutions is reported on Fig. |H Remark that 
for i?3 = 1, we have A(0) = 0; in fact we have A(p) = as a consequence of the fact 
that the positivity argument applies in the equation for the function A if Rq = R% . 

We manage to construct doubly-twisted solutions for generic values of the poten- 
tial, i.e. not NO-embedded solutions. From now on we will use the parametrisa- 
tion of the coupling constants formulated in terms of PH,Ph an d we will concentrate 
on the case = 1,Ph — 2, (3 — 7r/4; corresponding to the physical region for 
—3/2 < e 3 < 3/2. The profile of a doubly twisted solution is reported on Fig. [5j Let 
us stress that the functions A and a 3 are not proportional to each other for generic 
doubly twisted solutions. The non vanishing components of the electromagnetic fields 
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FIG. 3: The profile of a mixed solution corresponding to n = 2, m\ = 0, ni2 = 1. 
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FIG. 4: The energy, the charge Q an d parameters A(0), 03(0), /2(C)) are plotted as functions 
of the ratio R3 for Rq = 1 and (J3 = 0.017. 
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dA 
dp 



b 



dp 



5- : 



dA 



dp p 



(27) 
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These functions are represented on FigI5]for an embedded solution (solid lines) and for 
the doubly twisted solution of Fig. Note that the embedded solution has E p = 
since A(p) = 0. 



1.2 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

log(1+ p) 

FIG. 5: The profile of a generic solution corresponding to ph = 1, ph = 2, e = 0.1, (3 = 7r/4. 




0.01 0.1 1 10 100 

p 

FIG. 6: The profile of the electromagnetic fields associated with an embedded solution 
(solid) and a generic doubly twisted solution (dashed) are compared 
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Several features of the stationary twisted solutions are summarized on FigJT] and 
Fig. [HI On Fig. [7] some parameters characterising the evolution of a doubly-twisted 
solution (/2(C)) 7^ 0, #2(0) 7^ 0) into a simply- twisted solution (in the limit #2(0) — ■> 0) 
are presented. Although the effect on the energy is rather small, we have checked 
numerically that the energy slightly decreases while #2(0) increases. Finally, on Fig. 
IU the evolution of the limiting solution of FigJ7|for #2(0) = is reported for decreasing 
/ 2 (0). As expected, the energy of the twisted solution is lower than the one of the 
untwisted one. The critical value u b can be read on the picture. Although it would 
need more checks, these results suggest that doubly twisted string can have lower 
classical energy than the simply-twisted and the untwisted ones. A systematic check 
of this statement would need more numerical and/or analytical work. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



9 2 (°) 

FIG. 7: The evolution of some data for a doubly twisted solution with Rq = 2, R3 = —0.2 
and for pn = 2,Ph = 1> and e = 0.1, j3 = tt/4. 

Further aspects of the doubly twisted solutions is presented on Fig. O Here some 
features of the solutions are superposed for two values of the parameter R 3 with 
all other parameters fixed (namely Ro — 1, ujq — 0.01) and #2(0) labels the branch 
of solutions. The solid and dashed curves correspond to R3 = —0.5 and R3 = 1.5 
respectively. 

In the limit #2(0) — > the branch bifurcates into a simply twisted solution. In the 
limit #2(0) — >• 1 it ends up into a non localized solution. This figure also confirms the 
property that doubly twisted solutions have a lower energy than the simply twisted 
ones at least close to the bifurcating point where #2(0) = 0. In fact the decrease of the 
enegy is very small for the dashed curve but more significant for R3 < corresponding 
to the solid curve. After reaching a minimal value, the energy increases with #2(0). 

So far, all figures are presented with fixed value of e% = 0. However, as suggested 
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g 2 = , A = 
















K(0) 







0.0 0.2 0.4 0.6 0.8 1.0 



f 2 (0) 

FIG. 8: The evolution of some data for a simply twisted solution with Rq = 2, R3 = —0.2 
and for pn = 2,ph = 1, and €3 = 0.1, = n/4. 




in 12] , this coupling constant plays a role in the study of bifurcations. The analysis 
of the full bifurcation pattern for varying €3 (see Fig. 1 of 12]) in the present context 
is out of the scope of this paper due to the richness of the solution's manifold. Some 
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FIG. 10: The evolution of some parameters as functions of the coupling constant for 
R = 1 and #3 = -0.5 

information about the response of doubly-twisted solutions to €3 is presented on Fig. 
[TUl This graphic is obtained for R = 1,R 3 = —0.5, #2(0) = 0.5 and ph = l,p# = 2, 
(3 = 7r/4 (corresponding to the physical domain —3/2 < £3 < 3/2). For increasing €3, 
finite energy solutions stop to exist for €3 ~ 0.62. In this limit we observe that uo 2 = 0; 
as a consequence the various functions stop to be localized. For negative values of €3, 
the numerical analysis becomes very involved for 63 — > —3/2. Our results strongly 
suggest that, in this limit, the solution ends up into a simply twisted solution with 
respect to the doublet 0(2), i.e. it has /^(p) = and g%{p) 7^ 0. The fact that A(p) 7^ 
in this limit is related to the parametrisation (1141) adopted for this field in terms of 

LJ. 

E. Oscillating solution 

In j6| it was argued that timelike (or electric, i.e. with uj 2 < 0), finite energy 
solutions do not exist but no such solutions was reported. We failed to construct 
such solution in the case of one doublet but in the case of two doublets, the numerical 
techniques allowed us to construct such solutions. Refering to the asymptotic expan- 
sion of 0] it can be expected that several functions associated with timelike solutions 
oscillate asymptotically. We obtained indeed several such solutions. The profile for 
one of them is presented on FigJTTl The parameters are such that a 2 = cTq — o\ < 0, 
explaining in particular the oscillations of the function g 2 . Because this function be- 
haves asymptotically according to g 2 ~ sin(icxp)/ y/p, the term p(g' 2 ) 2 in the energy 
density leads to an infinite energy. 
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VI. CONCLUSIONS 

The lagrangian considered in this paper admits numerous types of physically rel- 



evant classical solutions. The most studied ones so far are the sphaleron [131 . [14 
constructed by using spherically symmetric solutions in the limit g' = and with an 
axial symmetry [2lJ for the physical value of the Weinberg angle. It is believed that 
they are related to baryon number violations processes which likely occurred copi- 



ously at some stages of the evolution of the Universe 15|, [la, [13. It has been known 



for some time that new solutions: the bisphaleron [18(, bifurcates from the standard 



sphaleron for a sufficiently large value of the Brout-Englert-Higgs-boson mass [18|, [19 



The recently discovered twisted semilocal strings jl, @] and the superconducting 
electroweak strings 22| open new ways of investigations of classical solutions in the 
MSSM. These solutions have applications in cosmology. The cylindrically symmet- 
ric ansatz used in this paper transforms the formidable system of Euler-Lagrange 
equations into a system of seven differential equations with boundary conditions and 
depending effectively on the five parameters characterizing the potential. The so- 
lutions are characterized by several integers labelling the winding numbers of the 
different components of the scalar fields relative to the axis of the symmetry. The 
solutions present a rich pattern of bifurcations occuring at specific values of the cou- 
pling constants which parametrize the underlying BEH-boson masses. Solutions with 
A = a 3 = f 2 = #2 = (NO-type) likely exist for all values of the coupling constants. 
However for sufficiently high values of the BEH-boson masses, solutions bifurcate 
from the NO-type solution and new solutions exist. They depend on time through 
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a non trivial phase factor and can be made z-depending (denoting z as the coordi- 
nate associated with the symmetry axes) by a Lorentz boost in this direction. Very 
reminiscent of the bisphaleron-sphaleron bifurcation, the static twisted strings have 
lower energy than the untwisted one. In contrast to the bisphaleron, which exist for 
Mn ~ 12Myy, twisted semilocal strings arise in regions of the parameter space com- 
pletely compatible with the experimental lower bounds of the BEH-boson masses. 
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